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. Abstract 

Most researches on fluid dynamics are mostly dedicated to obtain the so- 
lutions of Navier-Stokes equation which governs fluid flow with particular 
boundary conditions and approximations. We propose an alternative ap- 
Qh| proach to deal with fluid dynamics using the lagrangian. We attempt to 

develop a gauge invariant lagrangian which reconstructs the Navier-Stokes 
equation through the Euler-Lagrange equation. The lagrangian consists of 
! gauge boson field A^ with appropriate content describing the fluid dynamics, 

Q^ ' i.e. Aft = ($, —v). An example of applying the lagrangian to the interaction 

of fluid in a solitonic medium is also given. 
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1 Introduction 



The fluid dynamics still remains as an unsolved problem. Mathematically, a fluid 
flow is described by the Navier-Stokes (NS) equation [T]: 



where v is fluid velocity, P is pressure, p is density and p is the coefficient of viscosity. 

In principle, the study of fluid dynamics is focused on solving the Navier-Stokes 
equation with particular boundary conditions and / or some approximations depend 
on the phenomenon under consideration. Mathematically it has been known as the 
boundary value problem. The most difficult problem in fluid dynamics is turbulence 
phenomenon. In the turbulence regime, the solution for the Navier-Stoke equation 
has a lot of Fourier modes, such that the solution is untrackable numerically or 
analytically. It is predicted that the strong turbulence has 10 10 numerical operation 
j2]. This motivates us to look for another approach rather than the conventional 
ones. This paper treats the fluid dynamics differently than the conventional point 
of view as seen in some fluid dynamics textbooks. In this approach, the fluid is 
described as a field of fluid buch. We use the gauge field theory to construct a 
lagrangian describing fluid dynamics by borrowing the gauge principle. The Navier- 
Stoke equation can be obtained from this Lagrangian as its equation of motion 
through the Euler-Lagrange principles. 

2 Maxwell- like equation for ideal fluid 

The abelian gauge theory £7(1) is an electromagnetic theory that reproduces the 
Maxwell equation. To build a lagrangian that is similar with the abelian gauge 
theory, we should 'derive' the Maxwell-like equation from the Navier-Stokes equation 
[3] . The result can be used as a clue to construct a lagrangianan for fluid that satisfies 
gauge principle. Considering the Navier-Stokes equation Eq. for an ideal and 
incompressible fluid, 



dv 
dt 



+ (v ■ y)v = — vP - pV 2 v , 
P 



(1) 




(2) 



V -v = . 



(3) 



Using the identity v x (V x v) — V(|w 2 ) — (v • V)tf, it can be rewritten as, 




and then, 
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Putting the scalar potenstial $ = ^v 2 + — , the vorticity u = Vx? and the Lamb's 
vector I = uj x v, the equation becomes, 

dv _ 

— = -wxu-Vf 
dt 

= -I- V$ . (6) 
Imposing curl operation in Eq. (JHJ) we obtain the vorticity equation as follow, 

— = -v x (u x v) . (7) 



In order to get the Maxwell-like equation for an ideal fluid, let us take divergence 
operation for Eq. (J0J), that is 

d 

V • 1 = -V 2 $ = p . (8) 

Here we have used the incompressible condition, while by definition the divergence 
of vorticity is always zero, i.e. V • u = 0. Imposing again curl operation, we have, 

d _ _ - _ 

— (v x v) = -Vxl-Vx(V$), 

at 

v x r 

using the identity V x (V • <p) — 0. 



-Vxi, (9) 



Now, let us consider the definition of the Lamb's vector I = uo x v . Taking the 
derivative d/dt in the definition we obtain, 

dl duj _ _ dv . . 

d-t=m xv + u}X m- (10) 



Substituting Eq. ® and (0), we get, 



dl 



where, 



Vxw = aj+a-, (11) 



a = -, (12) 
j = _^v 2 $+ [V x (v-u)] v + uj x V($ + w 2 ) + 2 [(V x tf) ■ v] v. (13) 
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These results induce a series of equations, 



V-i = p, (14) 

v x r = -f , us) 

V-w = 0, (16) 

Vx^ = aj + a—, (17) 

that is clearly the Maxwell-like equation for fluids. If the fluid velocity is time 
independent, then I = — V$. This is the "electrostatic" condition. We use these 
results to develop gauge field theory approach for fluid dynamics in the next section. 



3 Bosonic lagrangian for fluid 

The correspondences of the electromagnetism and the ideal fluid can be written as 
follow, 

if <-> ut , 

«-> T , (18) 
A «-> v , 
<-> $ , 

where B is the magnetic field, E is the electric field, A is the electromagnetics 
vector, is a scalar function, Tj} is the fluid vorticity, I is the Lamb's vector, v 
is fluid velocity and $ is the scalar potential. The same as the electromagnetics 
field, we have a four vector = (0, A) which can be interpreted as the four vector 
for fluid dynamics, = ($,—£?). In the electromagnetics field, the scalar and 
vector potentials, <fi and A, are auxiliary fields. On the other hand, in the fluid 
dynamics the scalar potential $ = |-u 2 + V describes the kinetic energy of fluid, 
while the vector potential v is fluid velocity. Similar to the electromagnetics field, 
the lagrangian density has the form of [I], 

Ins = -\f^ + gj^, (19) 

where, 

•F mv = d^Ay d u A^ . (20) 

This Lagrangian obeys the gauge principles, i.e. it is invariant under a particular 
local gauge transformation, 

Afi > A'p = A^ + -d^O , (21) 

where 6 = 6(x) is an arbitrary real constant. It is easy to show that the lagrangian 
density in Eq. ()19j) is invariant under this transformation. 
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The equation of motion governed by this lagrangian can be derived using the 
Euler-lagrange equation in term of A^, 

&> 91ns - ®™ = o (22) 

After a straightforward calculation, we obtain, 

d u (d ll A u - d v d v A») -gj, = 0. (23) 

Now integrating it over x u and considering only the non-trivial relation as v ^ p 
gives, 

a A,-d,A = - 9 f ixo J^ a f ix ,J . (24) 
Since A4 = —v, A Q = $, d Q = d/dt and <9« = V. we have, 

-^-V$ = -sJ, (25) 



where jj = § dx Ji = — § dXiJ Q . Concerning the scalar potential given by $ = 
^v 2 + V, we obtain, 

d ~ * 1 — * 
- - -V |i/| 2 -VF = -gJ . (26) 

Borrowing the identity |v |w| 2 = (£T ■ V)tT + u x (V x v), we get, 

dv -, -> ~ t 

— + (v-y)v=-W-vx0-gJ, (27) 

where = V x v is the vorticity. This result reproduces the general NS equation 
with arbitrary conservative forces (VV). The potential V can be associated with 
some known forces, for example, P/p, (Gm)/r and r/(V ■ v). Here, P, p, G, v + rj 
denote pressure, density, gravitational constant and viscosity as well. 



4 Interaction between soliton and fluid 

In this section we describe an idea to apply the theory described in the preceeding 
section. We give an example on applying the theory to provide a consistent way for 
the interaction between soliton and fluid system. Soliton is a pulse-like nonlinear 
wave which forms a collision with similar pulse having unchanged shape and speed 
jO]. The wave equations that exhibit soliton are the KdV equation, the Nonlinear 
Schrodinger equation, the Sine-Gordon equation, Nonlinear Klein-Gordon equation, 
the Born-Infeld equation, the Burger equation and the Boussiness equation. Con- 
sidering the Nonlinear Klein-Gordon as follow: 
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The equation is a continuum version of that describes a propagation of molecular 
vibration (vibron) in a— helical protein jH]. The vibration excitation in the a— helix 
protein propagates from one group to the next because of the dipole-dipole interac- 
tion between the group. The wave is called the Davidov soliton 8J. Davydov has 
shown that in a— helical protein soliton can be formed by coupling the propagation 
of amide— J vibrations with longitudinal phonons along spines and that such en- 
tities are responsible for mechanism of energy transfer in biological system [8]. If 
a— helical protein immersed in Bio-fluid, then the phenomenon can be described by 
the interaction of soliton with fluid system. In standard technique in fluid dynam- 
ics, the problem will be done by solving of the Navier-Stokes equation and nonlinear 
Klein- Gordon simultaneously. 

In our current approach the problem is treated as follow. First, let us rewrite 
Eq. PHj) into four vector formalism, 

d^(j) - m 2 (j) + = . (29) 
Using the Euler-Lagrange equation, the lagrangian density is, 

I 771 ^ A 

I = - 2 {d,<P){d»<P) + —<f - . (30) 

In order to couple this lagrangian with the Navier-Stoke lagrangian in Eq. it 
is sufficient to replace the covariant derivative in Eq. (f3llj) [4 , 

= (d„ + igAJf . (31) 

The covariant derivative is invariant under local gauge transformation [6]. Then the 
interaction between soliton and fluid system obeys the lagrangian, 

11 m 2 A 

/ = — A F^ V + + ~Y<P 2 ~ 4,0 4 • (32) 



One interesting case is when we consider a static condition, i.e. d t f = with / is 
an arbitrary functions. Substituting A^ = ($, —v) into Eq. (|3*2"|) then the Lagrange 
density becomes, 

1 1 m 2 \ 

I = --(V x vf + -|(V - igv)<P\ 2 + - ^ • (33) 

The lagrangian is nothing else similar with the Ginzburg-Landau free energy 
lagrangian that is widely used in superconductor theory [7j. We have seen that the 
phenomenon of a— helical protein immersed in fluid similar with quantum electrody- 
namics for boson particle, while for static case it is similar with the Ginzburg-Landau 
model for superconductor. 
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Figure 1: Single soliton solution of the nonlinear Klein-Gordon equation. 



In order to perform an explicit calculation, suppose we have one-dimensional 
velocity in x direction v = (u(x),0,0) and <fi = <fi(x). Then the lagrangian in 
Eq.(jn31) reads, 

/ = leg - -g 2 u 2 <p 2 + - V 

2 x 2 2! 4! 

Substituting it into Euler-lagrangian equation we arrive at, 

d 2 (b , , , A 



(34) 



dx2 7Wn/ = 0, 



(35) 



where 7(x) 



m 



g 2 u(x) 2 . The equation is called the variable coefficient of non- 



linear Klein-Gordon equation. 

Further, we can consider a special case when the fluid velocity is constant, u(x) 
U, to obtain 



d 2 



dx' £ 



— 70 + a<p 







(36) 



with 7 = m 2 — g 2 U 2 and a = A/3!. To solve the equation, we can use a mathematical 
trick as follows. First multiply it by d(p/dx, 



d(f>d 2 (p d(p 3 d(p 
ax ax z ax ax 



(37) 



then integrating out over x and putting the integration constant as zero due to 
integrable condition lim^-boo = 0. Finally we obtain, 



dx 



a 

+ 2'- 







(38 



and it can be rewritten further as, 



a dx , (39) 



(f)(5 2 -0 2 )2 J V 2 

where 5 2 = 2j/a. Integration of the left hand side and solving the equation for 
provide the result, 



1 + e -Vf Sx e V f fa + e -Vf fa 
f^O = *sech ( yffe). 



coah(Vffa) V2 - (40) 

Thus, the solution for a homogeneous nonlinear Klein - Gordon equation is, 

<j){x) = Asech(Ax) , (41) 

where A = {12j)/X and A = (12 v / 37)/A 3/2 . This result is depicted in Fig. [TJ The 
figure shows that the soliton propagation will be damped by fluid. This theory also 
can be applied in turbulence phenomenon [3]. 



5 Conclusion 

We have shown an analogy between electromagnetics field and fluid dynamics using 
the Maxwell-like equation for an ideal fluid. The results provide a clue that we 
might be able to build a gauge invariant lagrangian density, the so-called Navier- 
Stokes lagrangian in term of scalar and vector potentials A^. Then the Navier-Stokes 
equation is obtained as its equation of motion through the Euler-lagrange principle. 
The application of the theory is wide, for instance the interaction between Davydov 
soliton with fluid system that can be described by the lagrangian density which is 
similar to quantum electrodynamics for boson particle. In the static condition, the 
lagrangian density is similar with the Ginzburg-Landau lagrangian. If the fluid flow 
is parallel with soliton propagation we also obtain the variable coefficient Nonlinear 
Klein-Gordon equation. Single soliton solution has been obtained in term of a 
second hyperbolic function. The result showed that the present fluid flow will give 
a damping in solitary wave propagation. 
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